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The statistical physics and dynamics of double supported bilayers are studied theoretically. The 
main goal in designing double supported lipid bilayers is to obtain model systems of biomembranes: 
the upper bilayer is meant to be almost freely floating, the substrate being screened by the lower 
bilayer. The fluctuation-induced repulsion between membranes and between the lower membrane 
and the wall are explicitly taken into account using a Gaussian variational approach. It is shown 
that the variational parameters, the "effective" adsorption strength and the average distance to 
the substrate, depend strongly on temperature and membrane elastic moduli, the bending rigidity 
and the microscopic surface tension, which is a signature of the crucial role played by membrane 
fluctuations. The range of stability of these supported membranes is studied, showing a complex 
dependence on bare adsorption strengths. In particular, the experimental conditions to have an 
upper membrane slightly perturbed by the lower one and still bound to the surface are found. 
Included in the theoretical calculation of the damping rates associated with membrane normal modes 
are hydrodynamic friction by the wall and hydrodynamic interactions between both membranes. 

PACS numbers: 



I. INTRODUCTION 

Double supported membranes are composed of two 
lipid bilayers superimposed on a solid substrate . The 
lower membrane (closer to the substrate) is weakly ab- 
sorbed, the upper one being bound via inter-membrane 
interactions. This system has a growing interest for 
physicists and biologists since it allows one to avoid, in 
part, the direct influence of the substrate on the upper 
membrane (pinning by defects, repulsion or direct attrac- 
tion by the wall) which becomes almost freely floating. 
Some recent experimental techniques allow the formation 
of such systems whether by Langmuir-Blodgett deposi- 
tion [IH3] or by deposing after rupture a giant vesicle on 
a single bilayer [H [5] . These double supported lipid bi- 
layers can then play the role of model membranes, the 
lipid and protein composition of which can be varied. 
Indeed, the extreme complexity of cell membranes mo- 
tivates the development of artificial ones [Zj that are 
more easily studied from a physical point of view. In 
the present case, the planar geometry facilitates the use 
of modern spectroscopy techniques to characterize these 
model membranes [U HI [5j [8] . Furthermore, two stacked 
membranes can model cell-cell junctions [9, 10 where the 
role of lipids and proteins can be investigated [TT] . 

Since liquid membranes such as lipid bilayers are highly 
fluctuating at room temperature (their bending modu- 
lus k is on the order of several k B T) [T3], a fluctuation- 
induced repulsion appears at finite temperature, leading 
to a destabilization of the system and eventually, an un- 
binding from the substrate at a critical temperature T u . 
The unbinding transition of membranes is a long stand- 
ing issue, first studied by Helfrich [T3], which has led to a 
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FIG. 1: Sketch of a double supported lipid bilayer. 



considerable amount of theoretical investigations on ho- 
mogeneous stacks of two or more membranes. Lipowsky 
describes, in a very nice review |14j . the various situa- 
tions where unbinding occurs, among which are the cases 
of only steric repulsive interactions with applied external 
pressure, or attractive interactions between membranes, 
depending on experimental conditions. 

These theoretical studies essentially focus on the order 
of the transition and the values of the critical exponents, 
mostly using Monte-Carlo numerical simulations [HI [16] , 
or group renormalization techniques [17H19) . Interest- 
ingly, an analytical solution can be obtained when con- 
sidering the unbinding of fluctuating strings (one dimen- 
sional objects) in two dimensions [50]. Indeed, this case 
is exactly mapped onto the derealization transition of 
a quantum particle in an external potential and the un- 
binding temperature T u can been computed exactly by 
solving the Schrodinger equation by transfer matrix tech- 
niques. 

Among these studies one can distinguish between sym- 
metric systems (two identical membranes or a bunch of 
identical membranes), and asymmetric systems where 
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the membranes are not identical-and eventually the ex- 
treme case of one membrane being a solid substrate (for- 
mally a membrane with infinite rigidity). Double sup- 
ported membranes belong to the latter category. The 
physics of asymmetric stacks studied by Monte-Carlo 
simulations [HI [21], or by analytical tools for one- 
dimensional strings |22) . turned out to be very rich, show- 
ing for instance, over a certain range of the parameters, a 
peeling process where successive unbinding temperatures 
appear, the upper membrane unbinding first. This has 
to be compared to the case of symmetric stacks where 
T u is independent of the number of stacked membranes. 
However, not much work has been done for real two- 
dimensional membranes embedded in a tridimensional 
space. 

Moreover, the major part of these studies consider the 
effect of bending rigidity in height-height fluctuations, 
without considering the microscopic surface tension <r, 
which might play an important role in supported bilay- 
ers. In such systems, a is related to the chemical poten- 
tial of the amphiphilic molecules (an energetic parameter 
associated with the full membrane area and not the pro- 
jected one) 23 . The value of this parameter is measured 
by fitting X-ray reflectivity spectra pQ, where it must be 
taken into account for wave-vectors q < Wcr/n, and in 
flickering experiments [23]. Moreover, according to the 
preparation method, the presence of pinning defects can 
induce lateral tension |10) . 

The surface tension always dominates bending rigid- 
ity at large distances, and thus drastically reduces the 
membrane roughness. It leads to a fluctuation-induced 
interaction which decays exponentially at large distances, 
instead of a power law in z~ 2 in the bending rigidity only 
case |25j . By comparing it to the direct interaction, such 
as van der Waals attraction and screened electrostatic 
attraction, we thus enter in the weak- or intermediate- 
fluctuation regimes. It has been shown, using renormal- 
ization techniques that, in theses regimes, the order of 
the transition can change and discontinuous transitions 
might occur [18 ] H9 ] [26] . 

The techniques presented above which include Monte 
Carlo simulations, renormalization group calculations, 
and numerical transfer matrix methods prove successful 
for calculating quantitative values for critical exponent, 
but they are not always easy to implement. Moreover, 
they often do not enhance our intuitive and qualitative 
understanding of the problem. In this paper, we develop 
a variational approach of asymmetric unbinding where 
the inter-membrane and substrate-membrane potentials 
are modeled by Morse potentials. The variational ap- 
proach provides an approximation that is analytically 
simple to implement and can provide a direct link be- 
tween quantitative calculations and qualitative pictures 
of the physics of the problem. This analytical work has 
thus similarities with the ones on single membrane un- 
binding, which use self-consistent [131 [27H29] and varia- 
tional [3U| calculations. Our goal is to understand the 
physical mechanisms and the role of the key physical 



parameters (adsorption potential, bending rigidity, sur- 
face tension) in the thermodynamics of supported stacked 
lipid membranes. What are the parameter ranges for 
which double supported bilayers are stable? In particu- 
lar, non-linearities of potentials are taken into account by 
the variational parameters, which are effective disjoining 
pressures and inter-membrane distances, and their de- 
pendence with temperature is explicitly studied. 

Fluctuating membranes are often studied by light- 
scattering experiments, which yield information about 
the bilayer dynamics and the associated damping 
rates |31j . Moreover, bio-membrane dynamics is cru- 
cial for the study of the diffusion of integral membrane 
proteins [32) . which is influenced both by fluctuations 
dynamics (projection of the motion onto a reference 
plane [33J [33]) and hydrodynamics [35 . The dynami- 
cal fluctuations of a single lipidic membrane in a vis- 
cous liquid had been studied in the early seventies by 
Kramer [36] , who showed that the damping rate is driven 
by the ratio between viscous damping and bending en- 
ergy. The presence of an external "obstacle" , such as a 
second membrane [HJ [37J or a solid substrate [35] , mod- 
ifies substantially the damping rate as soon as they are 
close enough to induce hydrodynamic interactions. Non- 
monotonic behaviour has been observed in the last case. 
In this work, we generalize these results to double sup- 
ported lipid bilayers and show how temperature and ad- 
sorption potentials modify the wave-vector dependence 
of damping rates. 

The first Section presents the variational approach in 
the simplest case of one supported bilayer and a com- 
parison is made with previous works. This variational 
approach is then applied to two stacked supported bi- 
layers, where phase diagrams and correlation function 
for height-height membrane fluctuations are computed. 
Finally, these variational parameters calculated in equi- 
librium serve to describe the normal modes and damping 
rates of the system. The flow velocity field is established 
for the first time in this geometry and hydrodynamic 
friction and hydrodynamic interactions between the two 
membranes turn out to be central to understanding this 
complex dynamics. A discussion of our results is given 
in conclusion. 



II. ONE MEMBRANE SUPPORTED ON A 
SOLID SUBSTRATE 

We first consider a single membrane lying on a hori- 
zontal solid planar substrate at position z — 0. Mem- 
brane local position vector is, in the Monge representa- 
tion, R = xx+yy + h(x, y)i where h(x, y) is the height of 
the membrane and rj — (h(x, y)} is the average distance 
between the substrate and the membrane. 

The Hamiltonian of a supported membrane is given 
by the usual Hclfrich Hamiltonian of a single fluctuating 
membrane [311 139j plus a term taking into account the 
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interaction between the membrane and the substrate: 
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where r = ix + yy is the projection of R in the plane 
parallel to the substrate, V = xd x +yd y and S is the pro- 
jected membrane area. The membrane bending rigidity 
is k (— 20 — bOksT for a lipid bilayer) and a is the mi- 
croscopic membrane surface tension. 

Several contributions arise in the potential V{z): at 
short distances it is dominated by excluded volume inter- 
actions and hydrophobic or hydration repulsion (entropic 
interaction associated to the presence of water molecules 
inserted between hydrophilic lipid heads) [27 . At inter- 
mediate distances, it is dominated by screened electro- 
static and van der Waals interactions which are essen- 
tially attractive. More precisely, it is usually admitted 
that, if the Debye screening length is small, the potential 
between the substrate and the membrane, assumed flat, 
is essentially the sum of two terms [T3] : a van der Waals 
attractive part between a flat membrane of thickness e 
and an semi-infinite medium 40 
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where Ar is the Hamaker constant between interfaces 
substrate/water and lipid/water (Ah — 1 — 10 x 10~ 21 J), 
and a decreasing exponential repulsive part which takes 
both steric and hydration contributions into account 
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where the values of the hydration pressure Phyd (1 — 10 x 
10~ 2 J.m -2 ) and Ah y d (0-1 — 0.3 nm) are not precisely 
known. 

In order to focus on the physical mechanisms of un- 
binding, we choose to mimic the adsorption potential by 
a simpler one, the Morse potential 
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which contains the two essential features described above 
(repulsion at short length scales and attraction at inter- 
mediate ones, see[2]). It is characterized by only three pa- 
rameters, its depth D, its range a, and the position of its 
minimum d. In principle, the following variational treat- 
ment could be done for the potential defined in Eqs. ([2|- 
([3]). However the Fourier transform of Eq. (J2|, needed 
in the following calculations, leads to singular parts that 
should be regularized using heuristic functions |30j . 

A priori, replacing the potential Eq. ([2| which varies 
as —z~ 3 for z 3> e ~ 3 nm by an exponential decreas- 
ing potential Eq. Q could modify the nature of the 
transition. Indeed, we move from the weak-fluctuation 
regime, |Vvdw| ^ \Va\, to the intermediate-fluctuation 
regime, | ^Morsel ~ |Vfl| (where Vg is the fluctuation- 
induced repulsion which decays exponentially at large 
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FIG. 2: Substrate-membrane potential as a function of the 
distance to the substrate: sum of van der Waals Eq. Q and 
hydration Eq. Q contributions (solid blue line) and Morse 
potential (dashed red line). Parameters for the Morse poten- 
tial are adjusted to have a good match at the minimum with 
the true potential (Phyd = 1.5 x 10 -2 J/m 2 , Ahyd = 0.1 nm, 
A H = 1.5 x 1CT 21 J, e = 3 nm, D = 1CT 5 J/m 2 , a = 2 nm -1 , 
and d = 1 nm). 



distances). However, it has been shown that in these two 
regimes, first-order transitions can occur [BSdS], an d m 
the intermediate-fluctuation regime their occurrence is 
related to the comparison of the two decay lengths (see 
below). 

At finite temperature, thermal membrane shape fluc- 
tuations induce an entropic confinement of the bending 
modes (which depends on n and a) and modify in turn 
the repulsive component of the interaction jT3l [27J [30] . 
Since the full calculation of the partition function Z = 
j D/ie~^[k] where /3 — (fc^T)" 1 , is untractable, we use 
in the following, a variational approach where the full 
Hamiltonian % is approximated by a trial Gaussian one 
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where the potential is harmonic with two variational pa- 
rameters, the spring constant A, and the shifted equilib- 
rium distance rj. 

The variational free energy reads in the Gibbs- 
Bogoliubov form 



-Fvar — Fq + [H — "Ho) 



(6) 



where /3Fo = — In Zq is the free energy associated with 
the variational Hamiltonian Q and the subscript refers 
to quantities calculated using Eq. {5}. The Gibbs in- 
equality ensures that F mr > -Foxact when F V&T is mini- 
mized with respect to the variational parameters. Their 
values will thus be determined by minimizing Eq. ^ with 
respect to A and rj. It must be emphasized, however, that 
we restrict our choice of variational Hamiltonians to the 
subclass of quadratic (and thus symmetrical) potentials, 
and by doing so, we obtain an approximate value of the 
true minimum of the free energy or equivalently the true 
values of the unbinding transition temperature and of the 
critical exponents. 

The calculations of Eq. ^ with the Morse potential, 
Eq. Q, have been done in Ref. |H] in the context of 
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DNA denaturation, i.e. for a one dimensional string in a 
space of two dimensions. In our case, the free energy Fq 
is 
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and the function g(x) = L dt/(x + 2t + t 2 ) is defined 

by El 
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In the case of vanishing surface tension cr — > 0, since 



l{x — > oo) ~ Eq. (|10|) reduces to the classical result 



of Hclfrich with no surface tension [T3J [3D] 
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Hclfrich closed the calculation self-consistently and found 
(it 2 )o = fi s d 2 /A with [i s 1/6 for a simple steric com- 
ponent at z = [13] . In our variational approach, the 
solution depends on the surface tension a and on the re- 
spective values of the Morse potential parameters, D, a, 
and d. 

Since F var = Fq+S(V)q — A ^j-, the minimization with 
respect to A yields A^^t = S 9 ^g2° as a l rea dy noticed 



by Podgornik and Parsegian 
implicit equation for A 
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Solving Eq. (13), which leads to the variational parame- 



ters A* and rj* , amounts to finding the lowest variational 
energy F*^ r = F var (A* , r?*) , solution of the problem. 

Let us denote by F vaI the variational free energy where 
rj is replaced by its expression in Eq. Q 
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where the function f(x) = f Q g(t)dt is 
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(15) 

The free energy of an unbound membrane which fluctu- 
ates freely in the bulk is given by Eq. (14 1 with A = 0. 
Thus, the membrane remains weakly adsorbed on the 
substrate as long as F var (A) < F var (0) = 0. We intro- 
duce the renormalized parameters 
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where the coupling parameter 3 is simply the rescaled 
second derivative of the Morse potential at the minimum, 
and 9 is a temperature scale. Since thermal fluctuations 
decrease the interaction between the membrane and the 
substrate, we expect that for a finite temperature, we 
have A(T) < 2a 2 D or A(T) < 3. 

From Eq. ([9|, one sees that the average height of the 
bilayer, 77, is proportional to the average mean square 
of height fluctuations (u 2 )q. Since g(x) is a monotonic 
decreasing function which tends to for large x and di- 
verges for x = 0, decreasing A, i.e. decreasing A or n 
or alternatively increasing cr, destabilizes the supported 
membrane. This increase of rj has been recently observed 
by heating double supported bilayers which makes n de- 
creasing (~ 200k B T in the gel phase to 1 — 3k B T at its 
minimum) [2J. 

It is interesting to note that naively one would ex- 
pect the membrane to desorb when the fluctuation con- 
tribution in Eq. ^ diverges, i.e. for A — > 0. Techni- 
cally speaking, this type of transition would be contin- 
uous. However, as said in the Introduction, discontin- 
uous transitions can occur. Indeed, clo se to the transi- 
tion we have I < 1 and Eqs. |9])-(10) yields rj — d ~ 
— 3/ (8ir)k B Ta/a ln(A/4). The variational free energy 
Eq. ( 14 ) can be written as 
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and thus comparing both exponentials and prefactors, we 
find an unbinding transition in two cases: i) a continuous 
transition for T/9 < 2 and 3 < 2T/9 < 3 C = 4, and ii) 
a discontinuous transition for T/9 > 2 and large 5 > 
2T/9 > 5 C = 4. 

More quantitatively, the unbinding tempe ratu re T u is 
defined by F var (A) = together with Eq. (13), which 
yields 
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FIG. 3: (a) Phase diagram of an adsorbed membrane in a 
Morse potential for a fixed value of the surface tension a. The 
solid (red) line corresponds to the rescaled unbinding temper- 
ature T u /0 vs. the coupling parameter H which measures the 
bare adsorption strength. It divides the phase diagram in 
two parts, where the membrane is bound (bottom) and un- 
bound (top). The dashed line correspon ds t o the system with 
vanishing surface tension given by Eq. (20 1, (b) Same phase 
diagram in the (H,T/T^) space, i.e. for constant adsorption 
strength, D, and varying a. The surface tension increases the 
unbinding temperature (dots) which behaves like l/\/H for 
< S < 4 (solid line) and reaches the asymptote 1 for large 



The solution of this system is found numerically and the 
phase diagram is shown in |3ja) keeping a fixed. At 
low temperature, the membrane is bound to the sub- 
strate and the binding region grows with the adsorption 
strength 3, which is in this case proportional to D. 

In the vanishing surface tension case, a — > 0, one finds, 
from Eqs. ( 18 )-( 19 1 and Eq. ( fl2| ), that the unbinding tem- 
perature is given by T u (a — > 0)/6> = 2/(7r v / e)- v / S, i.e. 



T* u =T u {a = Q) = 



2 VkD 
e akB 



(20) 



We obtain the correct scaling for T u 14J. Eq. ( p0|, in 
unrescaled units, corresponds to the dashed line in|3[a). 
Note that with the realistic parameter values given in [2] 
we find T u * =0.94 T room . 

The same phase diagram is re-plotted in[3]jb) using the 
parameters T/T* and 3. The adsorption strength D is 
kept fixed, and thus the parameter 5 now controls the 



surface tension. One observes that the unbinding tem- 
perature increases when 3 decreases, hence when the sur- 
face tension increases. Indeed, the system is entropically 
stabilized when a =/= 0, since more degrees of freedom 
are accessible to the bound membrane at a given tem- 
perature. At low 2, the unbinding temperature follows 
the limiting law T u /T* ~ Indeed Eq. (18) leads 

to T u — > 29 for A — > 0. This limiting form corresponds 
to the solid line in[3]jb). At large H, i.e. for a — > 0. 
T u /T* -> 1. 

The disjoining pressure, or mechanical stress |27j . is 
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S dd 3cv ' 

since dF vai /dd — d(VMoise(h))o/dd. Hence the varia- 
tional parameter A(T), plotted in [4] for various values of 
3, readily gives the disjoining pressure as a function of 
temperature, which can be measured experimentally. To 
compare this case to the vanishing surface tension case, 
we rescale the parameter A by the spring constant in the 
Morse potential 2a 2 D as A = A/(2a 2 D) = A/3. Equa- 
tion ( 19 ) becomes 
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Note however, that our variational approach, which be- 
longs to the "superposition of fluctuation-induced and 
direct forces" type of approaches, is not appropriate in 
this strong-fluctuating regime (with a = and thus 
|Vfl| 3> |^Morse|) 5 since, in this regime, only renormaliza- 
tion group approaches yields the correct (second) order 
of the transition [T4"] . 

and (b) are plotted A [Eq. (19)] and A 
as a function of T at constant a and con- 
stant D, respectively. First of all, one observes that A 
(and 77) depend on T which is a signature of the non- 
linear potential as for thermal expansion in solids. We 
recover the straight result that at T — >• 0, A(T = 0) = 
^Morse(^) = 2a 2 -D, i.e. the spring constant with no en- 
tropic induced repulsion Qb)]. Furthermore, at fixed T, 
A increases with 3 (or D) faster than 3, and A decreases 
slightly when 5 increases (or a decreases). This is an 
important result in the experimental context, since the 
surface tension of supported membrane may vary due to 
defects where the membrane is pinned or depending on 
the experimental protocol (e.g. washing processes or su- 
perficial pressure exerted at the edge). The variational 
mean height of the bilayer rj, which is the true order pa- 
rameter of the transition, is deduced following Eq. (13). 
On observes that, for 5 = 3, A = at the unbinding tem- 
perature (or 77 diverges). More generally, as suggested in 



Eq. (17), our approach yields a continuous transition (or 



second order) for s < 4 (solid line in the constant D 
phase diagram shown in^b)). For 5 > 4, A(T U ) is finite 
which makes the transition discontinuous (first order). 
The point (3 C = 4,T c /9 = 2) is thus a tricritical Lifshitz 
point, the surface tension a controlling the order of the 
transition. 



6 




FIG. 4: (a) Variation of the variational parameter A rescaled 
by the adsorption strength S as a function of rescaled temper- 
ature T/0 (a constant) for different values of 3 = 3, 30, 100. 
The parameter S is in this case proportional to the adsorption 
strength D. The fluctuation induced repulsion increases when 
T increases thus leading to a decrease of A, until we reach the 
critical point at T — T u . (b) Same as (a) with the temper- 
ature rescaled by the unbinding temperature at a — 0. 3 is 
thus proportional to a~ 2 and the case 3 — > oo corresponds to 
vanishing surface tension, Eq. (22 1. The associated average 



distance between the substrate and the membrane, rj, given 
by Eq. ( 13 1 is plotted in the inset. 



where ho = is the position of the substrate plane. The 
variational Hamiltonian reads 

1 f 2 

n = - / d 2 r]T{a(VMr)) 2 + K (V 2 Mr)) 2 

+A l [h l (r) - fn-i(r) - ^] 2 } 

= 1 f d 2 ru T (r)Uu(r) (24) 

2 Js 



with u T (r) = (hx(r) - i]i,h 2 (r) - r\i — r)%) and 



U= (-ctV 2 + kV 4 )I 



4 W / Ai + A 2 -A 2 
-A 2 A 2 



(25) 



where I2 is the identity matrix. In the eigenbasis the 
modes are decoupled with eigenvalues (of the last matrix 
in Eq. pi}) 
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and the Hamiltonian Eq. (24 1 becomes diagonal. The 
normal modes v T (r) = (i>i(r), v 2 (t)) are defined through 
u(r) = Pv(r) where P is a matrix of rotation of angle 
<t> E [0, f ] defined by 



tan <b = 



A 2 



A 2 - Al 
^2 2 



(27) 



One observes that \i > X 2 > where X 2 — for A% = 
or A 2 = 0. The case A\ = (<f> = 7r/4) corresponds to a 
free stack of two membranes remaining bound together, 
whereas A 2 = ((f) — 0) corresponds to the case of a 
second membrane free while the first one remains bound. 

By proceeding as above, the free energy Fq is simply 
given by the sum over the modes 



III. TWO SOLID-SUPPORTED MEMBRANES 

In this section we consider two membranes superim- 
posed on a substrate, with an adsorption potential as 
Eq. Q for the first membrane and a similar one for the 
inter-membrane potential. Since the range of the adsorp- 
tion potential is small (a" 1 ~ d, a few nm), the second 
stacked bilayer does not feel directly the substrate and 
is adsorbed only through the inter-membrane potential. 
We model these two potentials as Morse potentials Vi 
defined in Eq. Q with parameters Di, di and where 
z = l refers to the lower membrane and i = 2 to the 
upper one. The Hamiltonian of the system is 



+ y i [/i i (r)-/ li _ 1 (r)]} 



f 7(V^(r)) 2 + ^(V 2 ^(r)) S 
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where f(x) is defined in Eq. (15 1. Minimization of the 



variational free energy with respect to the four variational 
parameters Ai and r]i, yields the generalization of Eqs. ^ 
and (p| 



~cti{(ui - w i _ 1 ) 2 ) 



A l = 2a 2 D i exp [-a 2 ((u* - Uj_i) 



where 



2 k B T ^ ( 4Aj k 
27rer ^-^ 



,.9 



i=i 



V 2 



dX^ 
dA, 



(29) 
(30) 

(31) 



By replacing the r\i as a function of Ai and writing 
with a tilde dimensionless quantities, A4 = AAin/a 2 and 
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Aj = 4Ajk/(t 2 , the variational free energy F vax (Ai, A 2 ), 
corresponding to Eq. ( 14 ) for a single supported mem- 
brane, is 



167TK 



/(A X ) + /(A 2 )- 



OA 



exp 



27TCT 



(32) 



where Ai, 2 (Ai, A 2 ) are given in Eq. (26 1. 

The unbinding of the adsorbed stack can a priori fol- 
low three different scenarii : i) the upper membrane des- 
orbs and diffuses freely in the bulk while the lower one re- 
mains adsorbed (A 2 — 0); ii) the stack as a whole evolves 
freely in the bulk (A\ — 0); or hi) the two membranes 
are completely de-stacked and desorbed. Hence, to study 
quantitatively this unbinding transition, the variational 



free energy, Eq. ( 32 1 , should be compared to the varia- 



tional free energies of the system in the three configura- 
tions described above 



AF vaI (A 1 ,A 2 ) = F var (A u A 2 ) - F v: 



(33) 



where -FVar is the variational free energy of the (possi- 
bly partially) unbound system. Since the variational free 
energy of a freely fluctuating membrane in the bulk has 
been chosen as the reference of energy, these free ener- 



gies, 



are res pec tively i) Eq. ( 14 ) at its minimum, 
F var (A*)} ii) Eq. |32]) with Ai = 2A 2 and A 2 = at its 
minimum A% 




FIG. 5: Phase diagram of two adsorbed membranes in a Morse 
potential, (a) Hi = 30; (b) Hi = 100. The colored dots cor- 
responds to the rescaled unbinding temperature T u /6 of the 
adsorbed stack vs. the coupling parameter H2 which mea- 
sures the strength of the inter-membrane potential. Solid 
lines (dashed for unstable case) correspond to the unbind- 
ing of the single membrane when the upper is unbound (red) 
and of the two membranes in the free stack (blue). The dia- 
gram is thus divided in four regions, where the stack is bound 
(region 1), the upper membrane desorbs (2), both membranes 
are unbound (3), and the stack is unbound (4). 



^var(^) 



167TK 

-/3D 2 exp 



f(2A* 2 )~2A* 2 g(2A* 2 



(34) 



and hi) F var (0) = 0. 

H ence , the unbinding temperature, T Ul is solution of 
Eq. (|30|) and AF var (Ai,A 2 ) = 0, i.e. the values of A x , 



A 2 and T u are solutions of the three equations: 



E 



/(AO -£ 5(^1^-^ 
.. \ a A. J u 

3 = 1 1 



f{At)-g{A\)A\-4- u A\ 
f{2A* 2 ) - 2g(2A* 2 )A* 2 " 




v A* 

T„ 2 



(35) 



Ai = ^exp 



i€{l,2}.(36) 



The phase diagram is shown in [5]ja) for Si = 30 and in 
[5jb) for Si = 100 (ai — a 2 ), and S2/S1 varying between 
and 10. The surface tension a and bending modulus k 
of both membranes are kept fixed, which sets the parame- 
ters Si to be proportional to the two adsorption strengths 



Di. Four distinct regions appear, related to the four free 
energies defined above. 

For low S 2 < 2Si and by increasing temperature, the 
stack is progressively peeled up, the upper bilayer un- 
binding at a temperature lower than for the bilayer close 
to the substrate (region 2). However, for intermediate 
values of S 2 [2Si < S 2 < 235.5 in 5^a) and 555 in[5jb)], 
the stack unbinds completely at the transition, defining 
one unique unbinding temperature T„(S 2 ). One then en- 
ters in region 3 where both membranes are unbound. 
This unbinding transition temperature of the adsorbed 
stack is larger than for a single adsorbed membrane. In- 
deed, the upper membrane reinforces the adsorption of 
the lower one by "squeezing" it to the substrate. More- 
over, the region in the phase diagram where the stack 
exists is larger than when there is no adsorbing surface. 
This is due to the slight decrease of the first membrane 
height fluctuations induced by the presence of the wall. 
Finally, for large enough S 2 values, a fourth region ap- 
pears (region 4) where the stack unbinds as a whole, 
the two membranes remaining bound together because 
of their strong mutual attraction. 

In [6] are plotted the variational parameters as a func- 
tion of T/8, for three different values of S 2 /Si = 1,5 and 



H 2 =30 H 2 =150 S 2 =300 




6 I g 1 6r 1 6 

a(ri-cO 5 5 5 




0.5 1 1.5 2 2.5 3 3.5 1 2 3 4 1 2 3 4 5 6 



779 T/Q 7/9 

FIG. 6: Top: Variational parameter A/3 = A/(2a 2 D) vs temperature T/0 for S2 = 30, 150, 300 (Hi = 30). Red dots correspond 
to the lower bilayer, blue ones to the upper one. For comparison sake, solid lines correspond to one supported membrane for 
S = 30 and a free stack of two bilayers for H2 = 300. Bottom: Corresponding average distance between membranes, 0(77 — d), 
given by Eq. Q, vs. T/6. 



10 (Si = 30), corresponding to the three regimes just de- 
scribed and shown in [5^ a) . Clearly, one observes that 
for the upper membrane, A%(T) increases when S2/S1 
increases, as expected. For S2/S1 = 1, A\{T) > A 2 (T) 
which is signature of the increase of height fluctuations: 
the upper membrane fluctuates more than the lower 
thanks to the fluctuations of the lower membranes which 
add up with its own ones. For 3 2 /Si = 5 we have 
Ai(T) ~ A 2 (T), which are both larger than above, and 
the associated interlayer distances are almost the same. 
This is the reason why the stack unbinds completely at 
the transition. Finally, for S 2 /Ei = 10, Ai(T) < A 2 (T) 
and the transition of the stack is reached when Ai(T u ) = 
0. We see that the adsorption strength of the substrate 
Si compared to the interlayer potential tuned by S 2 is 
central in computing membrane fluctuations. 

More interestingly is the nature of the transition as 
a function of S 2 : for low S 2 < 2Si and very large S 2 
the transition is continuous occurring respectively for 
A 2 (T U ) = and Ai(T u ) = 0, whereas for intermediate 
values of S 2 corresponding to the full unbinding of the 
stack (from region 1 to region 3) the transition occurs 
for finite values of Ai, i.e. it is discontinuous. This can 
be related to the unbinding of one bilayer with a direct 
interaction including a potential barrier [141 119] which 
exhibits discontinuous (or first order) unbinding transi- 
tions. In our case, the upper membrane induces such a 
potential barrier felt by the "squeezed" membrane. 

Once the variational parameters, Ai, are determined, 
one has access to the fluctuations of the two membranes 
and the height-height correlation functions. The varia- 
tional Hamiltonian being Gaussian, the structure factor 



is 



(^(q)^(q')) = (27r) 2 ( 5 i ^(q + q') 



k B T 
Aj + aq 2 + nq 4 



(37) 



and the three height-height correlation functions 
Cjj(r) = (iij(r)uj(O)) are given in Appendix |A| In par- 
ticular, we find 







cos 2 (j> < 







k B T 
2na 



I \ kBT J. ■ 

\U1U2) = cos sin < 

2ira 



k B T 
2na 



2 , 4AiK 

sin <p g 



+ sin 4> 9 



4A 2 k 

T 2 " 



cos <f> g 




which are plotted 



Ml) 



-<t)\ 



Note that 



see Eq. (31 



in|6] For S 2 = ai the fluctuations o 



3ai 



and 



] are already plotted 
the upper membrane 
increase from one to two supported bilayers. In this case, 
fluctuations of the lower membrane add up to the upper- 
membrane ones leading to a substantial increase for T > 
29. Membrane 1 is thus very slightly perturbed by the 
presence of membrane 2 (compare lines to dots in [6] and 
[7| . Finally, the correlations between both membranes are 
very low. When S 2 increases, correlations between both 
membranes increase and their correlation functions are 
almost identical (for S 2 /Si = 5 and 10). 

In principle, the same calculations can be done for n 
membranes (see Appendix [B| with eigenvalues Ai where 
i G {0, ...,n — 1}. The variational equations are the 
same as Eqs. p5|)-(36), where the rhs. of Eq. (35) now 
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FIG. 7: Height-height correlation functions given by Eqs. ( 38 1 , (391, and (401 (a = oti = 012) as a function of temperature for 
increasing H2 (same values as in|6|. Solid lines correspond to one supported membrane (H = 30) or a free stack of two bilayers 



(S 2 = 300). 



contains all the free energies for free bundles made of j 
membranes (A n _j±i = 0) and adsorbed bundles made 
of n — j membranes. Eigenvalues and eigenmodes in the 
simplest case where all the variational parameters are 
equal, are given in Appendix [B] 

In this section, we have implicitly assumed 012 = oti. 
However, one might expect to have actually different po- 
tential ranges ai > ai since the inter-membrane poten- 
tial is at large distances in z~ 4 instead of z~ 3 [13]. It will 
introduce two temperatures 9\ > 82 and the results will 
remain qualitatively the same. 



IV. DYNAMICS OF TWO SUPPORTED LIPID 
BILAYERS 

Dynamics of supported membranes are governed by 
Langevin equations for height displacements of mem- 
branes 1 and 2, hi(r) and /12(f), written in q-space 



^^=- A - (q) ^(= q ^ +Cl(q ' <) 



(41) 



where Q (<L t) is the Fourier transform of the ran- 
dom noise obeying the fluctuation-dissipation theorem: 
(Ci(q,*)C*(q',*)) = 2(27r 2 )k B TA tJ (q)S(t - f )*(q + q'). 
The damping matrix A(q) takes hydrodynamic interac- 
tions into account, both along the membrane and be- 
tween adjacent surfaces (substrate-membrane 1 or mem- 
brane 1-membrane 2). 

It is known that for a single membrane in an infinite liq- 
uid, the mobility is simply A(g) = (4/ig) _1 and Eq. (41) 
leads directly to the damping rate [3T] 



7o(g) = (nq 3 + aq)/4/2. 



(42) 



It is the ratio of the energy driving fluctuations and the 
viscous damping. However, in the case considered in [T] 
we expect that both the presence of the second mem- 
brane and the solid wall, which imposes a no-slip condi- 
tion, will substantially modify Eq. (42 1. In the following, 
we compute A(q) in a simpler geometry, also shown in[l] 
(dashed lines), where the membranes are supposed pla- 
nar, located at z = rji and z = 772. 



A. Hydrodynamics of two supported membranes 

The velocity flow field v(R) and pressure field p(R) 
in the bulk, are found using the Stokes equation for an 
incompressible fluid of viscosity fi 



V • v 

/A7 2 v - Vp 



(43) 
(44) 



with the following conditions at the boundaries 



p(r, z — > 00) 


= Pq 


Vr 


(45) 


v(r, z — > 00) 


= 


Vr 


(46) 


v(r,0) 


= 


Vr 


(47) 


v ( r ^r 2 ) 


= v(r 




(48) 


V|| • vn(r, 771,2) 


= 




(49) 



Eq. (47) imposes the no-slip condition at the substrate, 
and Eq. ( 49 ) ensures the membrane incompressibility. 



The two fluctuating membranes impose normal forces 
/1 2 at z = 771,2 which are balanced by fluid stress jumps 



fi,2 = -<7zz(r,Vi,2) + a *z(r,r)i,2) 
where the stress tensor is cr i? = —pSij + /i(dji>i 



(50) 



From Eqs. ( 43 )-( 49 ) , one finds easily that 
A,2 = <^(r, 771,2). 



(51) 



Due to the linearity of Eqs. ( 43 )- ( 44 ) , one expects a linear 
relation 



(52) 



between normal forces applied on membranes /1.2 and 
the z-component of flow velocity at the membranes, 7>i,2j 
which defines the resistance matrix 4/igL. Following 
Brochard and Lennon [53] and Seifert [3S], one seeks for 
solutions of the type v\\ = v x {x 1 z) — j '{z)e lqx ~ lt , and 
v z (x,z) = g{z)e l( i x -^. From Eqs. |43}-((47} we find, by 
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writing z = qz and fji = qrji 
gi(z) = A\ sinh z — z cosh 



fji sinh f\\z sinh z 



pi(z) = 2/j.qAi 



cosh ; 



(54) 



sinh fji + fji cosh fji 
( fji sinh fji sinh z, 
\ sinh fji + 771 cosh fji 

g 2 {z) = {A 2 +B 2 z)e i + {C 2 + D 2 z)e- i (55) 

P2 (S) = 2fiq(B 2 e s + D 2 e- i ) (56) 

53(5) = (C 3 +^)e- 5 (57) 

p 3 («) = 2 M D 3 e- s (58) 



Finally by solving the system of four equations Eqs. |48|)- 
(49 1, we determine the four coefficients A 2 ,B 2 , C 2 and D 2 



and inserting the result in Eq. (51 1, we find 



where 6 — q8 = q(r/ 2 — 971) ((5 
distance between membranes), 



L(fji,5) = A{fj 1 )+B(5) (59) 
? ?2 — Vi is the average 



A(r?i 



r/i +cosh r/i sinh 771 
sinh 2 771 — r/ 2 





(60) 



and 




(5+cosh 8 sinh <5 
sinh 2 S—S 2 



sinh 2 8 — S 2 



sinh 8+8 cosh (5 

sinh 2 8— 8 2 
1 £+cosh 6 sinh <5 
h sinh 2 S-8 2 



(61) 



The damping matrix in Eq. (41) is nothing but the in- 
verse of the resistance matrix 

A(q,r,i,r, 2 ) = ^L- 1 (qr ) i,q5). (62) 

The resistance matrix 4:fj,qA(fji) corresponds to the fluid 
friction induced by the substrate where the no-slip 
boundary condition applies, whereas the symmetric re- 
sistance matrix 4/j.qH(5) corresponds to the mutual hy- 
drodynamic interactions between both membranes. 

In the limit 771 — > ±00, the fluid friction disappears 
(A — > 0), and we are left with the symmetric resistance 
matrix 4/igB(<5), the large and low q limits of which were 
studied by Brochard and Lennon [531 ISZ] • On the other 
hand, in the limit 5 —¥ 00, membranes decouple since 
B —> I2 and we find the result of Seifert for a single 
membrane close to a substrate [38] . This is also the case 
in the limits fji — > or S — > where the matrix L becomes 
a scalar 



L{x) 



+ cosh x sinh x 
sinh 2 x — x 2 



(63) 



where x — 5 or x — fji respectively. The linearized lu- 
brication approximation x <C 1 in Eq. (63 1, i.e. L(x) = 



3/x , corresponds to a flow parallel to the substrate 



B. Damping rates 



Solving the coupled Langevin equations Eq. (41) 



amounts to finding the eigenvalues and the eigenmodes 
of the damping matrix 7(g) = \(q)\J(q) where 



U(g) = 



_ / aq 2 + nq 4 + A x + A 2 -A 2 



-An 



aq 2 + nq 4 + A 2 



(64) 



By writing the two eigenvalues 71(9) and 72(9), Eq. (41 ) 
becomes, in the eigenbasis 

= - 74 (gH(g,t) + (Rr%-Ci(<7,*) (65) 

where R is the transformation matrix, u = Rw. The 
solution for the temporal correlation function of Wi(q,t) 
is 

(w t (q, t)wi(-q, 0)) = ( Wl (q, 0) a )e-*»(«)* (66) 

and the correlation function (v,i(q, t)ui(— g, 0)) can then 
be computed coming back to the initial coordinates and 
using 



(Mg)| 2 ) = 



sin 



cos" 



Ai + aq 2 + nq 4 X 2 + aq 2 + nq 4 
1 



(67) 



(u 2 (q)ui(~q)) = sine/) cos <j> 



1 



A 2 + aq 2 



Kq 1 



Ai + aq 2 



nq" 1 



(68) 



The eigenvalues 7j(g) are increasing functions (see [8]) 
where limiting forms for q — > show a quadratic be- 
haviour 



7i(g) 



A 2$ 3 2 , + Vl 2 

-q and 72(g) = -^ttt - g 



12/i 



I2fi 



(69) 



For q — > 00 membranes decouple and we simply recover 
the free membrane damping rate Eq. ( 42 ) 



71(g) = 72(g) = ^g 3 - 



(70) 



Note that without taking into account hydrodynamics, 
the damping rates 7i(g) are 



woffle K 1 + ac T + A, 



7i 



l (g) = 



4yug 



(71) 



where A^ are the eigenvalues given in Eq. (|26j). Hence 
hydrodynamics induce a large decrease of the damping 
rates for small and intermediate values of q, which is 
expected due to the no-slip condition for the flow velocity 
at the substrate that slows down hydrodynamics. 
Introducing the elastic decay length 



(72) 
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10 7 J.m 2 which are reasonable values [TJ[TD]. Note that 
close to the unbinding transition for T ^ T u , the average 



FIG. 8: Log-log plot of the (dimensionless) damping rates 7; 
versus adimensional wave-vector q£ for Hi = 30 and H2 = 3 
at temperature T = 1.5 # (di = 2 nm and d,2 = 3 nm). Thick 
solid lines correspond to the full damping rates 74 (i = 1 in 
red and 2 in blue) whereas thin solid lines are damping rates 
without hydrodynamics Eq. (71 1, and the dashed line is the 
asymptotic result (g£) 3 - 



allows us to define the relevant dimensionless parameters 
in the dimensionless damping rates 7^ defined as 



1% <£,,M,M, 



7i(<?) 



(73) 



where the A\ and r\i are fixed for a given couple of 
(Sx.32). 

Damping rates Eq. (73 1 are plotted in [8] for (Si,S 2 ) = 
(30,3) at relative temperature T/6 = 1.5. Following 
the previous Section, it corresponds to variational pa- 
rameters equal to A\ = 18.6 and A 2 = 0.4. From the 
phase diagram in [5j we are thus in region (1) close to 
the unbinding of membrane 2 [a(r)i — d±) — 0.72 and 
01(772 — d 2 ) — 4.76]. However by choosing d\ = 2 nm 
and d 2 = 3 nm, this yields ?7i/£ = 0.07 and S/£ = 0.3, 
and the substrate-membrane 1 and inter-membrane dis- 
tances are small enough to induce a sensible slowing down 
at small and intermediate wave-vectors q. Indeed, when 
hydrodynamic interactions are neglected (thin lines in[8| 
damping rates are larger by 4 to 6 orders of magnitude: 
for q — > 0, 7j voHI (?) reduces to Aj/4/xg which should be 
compared to actual values in Eq. ( [69] ) . Moreover, the 
damping rates of the two normal modes are less different 
for low q than expected without hydrodynamics. From 
Eq. ( 69 ) the low q values of the damping rates 7$ (<?) are 
controlled by Aiijf for mode 1 and A 2 S 3 for mode 2. Usu- 
ally the values of these quantities are quite close, which 
can be explained by the fact that an increase of Ai im- 
plies a decrease of 77, as shown in the previous Section 
[see for instance Eqs. ( 30 )-( 31 1] , and the decoupling of 
the two modes occurs only at intermediate q (see [8j) . 

At room temperature, the values used in [H] yields £ = 
56 nm (for a = 7 x 10" 5 J.m" 2 and n = 30fc B T) with 
a7 l = aW 1 = 2.5 nm, D 1 ~ 10" 6 J.m" 2 , and D 2 ~ 



square value of height functions [Eq. (40])] takes large 
values, the approximation of almost planar membrane 
fails, and dynamical renormalization techniques should 
be applied. 



V. DISCUSSION 

The variational approach that we have developed to 
describe the statistical physics of double supported mem- 
branes is a complementary approach to numerical Monte- 
Carlo simulations and numerical transfer matrix meth- 
ods. It allows an analytical determination of the un- 
binding temperatures, and the variational parameters, 
the effective spring constants Ai(T) and the substrate- 
membrane and inter-membrane distances r)i(T), at any 
temperature and any values of the dimensionless adsorp- 
tion strengths (Si, 22). 

This is central for analyzing X-ray specular reflectivity 
spectra since the Ai(T) are necessary to fit the experi- 
mental structure factors. In the experiments by Daillant 
et al. pQ , the measured effective spring constant of the up- 
per bilayer A2, which has a large effect on the spectrum 
shape, was found surprisingly weak, A 2 ~ 5 x 10 10 J.m -4 , 
much lower than the bare spring constant of the potential 
estimated at 2a 2 D 2 ~ 2 x 10 14 J.m -4 . With these values 
we find A2/E2 — 2 x 10~ 4 , a value extremely low suggest- 
ing that the system is closed to the unbinding transition 
for the upper membrane. Inserting their fitted parameter 
values for the surface tension and the bending rigidity 
(cr = 5 x 10~ 4 J.m~ 2 and n ~ 1.5 x 10 -19 J), we find 
S 2 ~ 200. Assuming for the first membrane Si ~ 100, 
the phase diagram [5jb) shows that this would lead to 
T ~ 56*. With a ~ 2 nm" 1 we thus find T ~ 300 K which 
is roughly the temperature of the experiment. Although 
this is a crude estimation, we find the correct values, 
even if the parameters A^ depend in a complex manner 
on both inter-membrane and substrate-membrane poten- 
tials. Note that more generally, even if the transition is 
not perfectly described within our variational approach, 
the error made on the value of the unbinding temperature 
T u is small since Ai(T) varies rapidly close to T u . 

In recent experiments by Stidder et al. [13] , unbinding 
of the upper membrane in DPPE double stacked bilayers 
with 10 mol% of cholesterol has been observed between 
52.2 °C and 56.4 °C by neutron reflectivity. Their results 
suggest i) a discontinuous unbinding transition since they 
observe both an hysteresis and a larger but finite rough- 
ness, \J (u 2 ) = 0.6 ±1.2 nm, at the transition, and ii) 
that the presence of 10 mol% of cholesterol modifies the 
direct interaction potential and decreases the unbinding 
temperature. We qualitatively obtain such unbinding 
transitions for S2 J$ 2Si but values of the surface ten- 
sion and the direct adsorption potential are needed, to 
do a quantitative comparison. Nevertheless, these exper- 
iments allow one to hope that a detailed characterization 
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of double lipid bilayers will be achieved experimentally 
in the very near future. 

When compared to numerical results of Refs. [T51 I5T] 
where only the case Si = S 2 was explored, we also find 
sequential thermal unbinding transitions: for Si = S2 = 
30 the upper bilayer unbinds at T u (2)/6 = 2.6 and the 
remaining bilayer then unbinds at T u {\)/9 = 3.3 (see [6] 
and [4|. Moreover, membrane fluctuations increase when 
moving away from the substrate, 772 — e?2 > 171 — d\. We 
show that this behaviour is valid whenever S2 < 2Si. 
However, for higher values of S 2 , one observe the re- 
verse case where the stack of two membranes unbinds 
as a whole. This has been shown for bundles of strings 
in two dimensions [22 . Our phase diagram ^ resem- 
bles the one of Ref. [55] for strings. A way to compare 
them is to change the coordinates of the phase diagram 
into Si i2 (T/6>) 2 oc T 2 /(kDi i2 ) in order to eliminate the 
surface tension (which does not exist in the string prob- 
lem). Note that, for strings, the parameter k becomes an 
elastic (stretching) parameter. The qualitative diagram 
is then similar with a region corresponding to the peel- 
ing process and another one where the bundle desorbs as 
a whole. However, contrary to these studies, we find a 
range of parameters S2 for which the double supported 
membranes de-stack completely at the transition, this 
region being larger and larger when S 2 increases. The 
richness of the phase diagram is related to the fact that, 
contrary to Refs. [13 [5TJ [55] we consider bilayers with 
a microscopic surface tension, which is the more general 
case. 

In this work, we assumed the substrate to be flat and 
did not considered the substrate roughness which leads 
to pinned or hovering binding states for the first mem- 
brane [UJ . How does the influence of the quenched sub- 
strate disorder propagate to the second membrane is an 
important issue. One might suggest that this effect for 
the upper membrane is somehow annealed by thermal 
fluctuations of the lower one, and the pinning effect is 
then less pronounced. Anyway, the influence of substrate 
roughness is central from an experimental point of view 
and deserves a quantitative study. 

From an experimental perspective, we now discuss the 
biophysical interest of designing supported double (or 
even multiple) lipid bilayers. Ideally, for the upper bi- 
layer to deserve to be qualified as "nearly floating" and 
weakly perturbed by the substrate, both 6 and 1 /A2 must 
be large so that this bilayer fluctuates nearly freely, while 
still being conveniently observable by spectroscopy tech- 
niques. One can in principle obtain a single adsorbed 
bilayer with the same properties, by approaching its un- 
binding transition. However, when the experimental goal 
is to model a biological situation (e.g. plasma membrane 
or cell-cell junction), both membrane composition and 
temperature are imposed by the biological context. This 
reduces the possibility to design a single "nearly float- 
ing" bilayer. By contrast, our study suggests that such a 
regime can be reached by stacking two (or more) bilayers, 
because thermal fluctuations are enhanced in the upper 



one. Our goal is to help to anticipate in which regime of 
parameters this situation is likely to occur. 

Another experimental interest of studying thermal 
fluctuations, both at the equilibrium and dynamical lev- 
els, lies in the possibility to study molecular diffusion in 
supported bilayers, for example by single molecule track- 
ing. On the one hand, fluctuations affect apparent dif- 
fusion coefficients because the projected area is different 
from the real membrane one; the relationship between 
apparent and real diffusion coefficients depends on the 
relative timescales of membrane fluctuations and diffu- 
sion [35]. On the other hand, it has been recently shown 
by numerical simulations |35j that membrane elasticity 
decrease protein mobility, essentially due to the viscous 
drag of the deformed membrane around the protein, but 
that membrane fluctuations slightly increase protein mo- 
bility 33J. However such effect is weak compared to vis- 
cous losses [55] , It will be useful to quantify these effects, 
related to the amplitude of fluctuations, in the case of 
stacked bilayers. This work is currently in progress [38] , 
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Appendix A: Correlation functions for double 
supported bilayer 



By performing an inverse Fourier transform of Eq. ( 37 ) 



and coming back to initial coordinates, one finds the fol- 
lowing height-height correlation functions 



C n (r) = ( Ul (r)«i(0)) 



k B T 
2ir 
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cos 



dq 



qM<i\ r \) 



dq 



Ai + aq 2 4 
A 2 + aq 2 + nq 4 



(Al) 



C 12 (r) = (ui(i> 2 (0)) 
ksT cos </>sin< 



+ 



2tt 
dq- 



dq 



qMq\A) 



Ai + aq 2 



qJo(q\r\) 



10 A 2 +cr<2 2 
C 22 (r) = (u 2 (i> 2 (0)) 
k B T 



nq H 



Kq* 
(A2) 



2n 



- COS 



sin 



dq 



qMq\ t \) 



dq 



Ai + aq 2 
gJo(g|r|) 



Kq 



A 2 + aq 2 



Kq* 



(A3) 



where Jo is the Bessel function of the first kind, and both 
<t> and Ai depend on the variational parameters Ai(T) and 
A 2 {T), following Eqs. ([26} and ((27). 
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Appendix B: Stack of n supported membranes 

In the case of n stacked membranes on a substrate 
at z = 0, the full Hamiltonian can be written as the 
tensor product of the usual Helfrich Hamiltonian of a 
single fluctuating bilayer membrane and the Hamiltonian 
which describes interactions between membranes: 

phonon H~l Helfrich onon[V 

(Bl) 

where "Hncifrich is the classical Helfrich Hamiltonian of a 
free membrane, Eq. ([I]) with V = 0. The phonon Hamil- 
tonian is described by a succession of harmonic springs 
in the z-direction 



1 

3=0 



(B2) 



Taking thermal expansion into account would require a 
variational approach where the At arc different and com- 
puted variationally from, for instance, Morse potentials 
as defined in Eq. Q. However, even if possible in princi- 
ple, this approach is hardly tractable in practice. In the 
following of this appendix we consider all the Ai equal, 
Ai=A,Vi. 

Since the boundary conditions for the phonon modes 
are: hp(r ) — and the last membrane free, we find from 
Eq. (B2) the following eigenmodes and eigenvalues 



m = - 



4 sin 



2 . / 21 ■ 

. sin 

V2n + 1 \2n 

2 (21 + 1 tt" 



1 



IT J 



V 2n + 1 2 



(B3) 
(B4) 



with < £ < n - 1. For n = 2 we find A = 4sin 2 (7r/10) 
and Ai = 4sin 2 (37r/10) which corresponds to A 2 /A 2 and 
\i/Ai respectively, in Eq. d26j) (with A 1 =A 2 =A). 



The Hamiltonian (Bl I is diagonal in the Fourier (£, q) 
space, where £ is the index for the deformation modes 



in the z-direction and q is the planar wave-vector in the 
(x, y) plane: 



H[h(t 



i r J 2 n—1 

' q)] = 2 J |^£(^WW)l^,q)| 2 

(B5) 

where we define the Fourier transform of height fluctua- 
tions hj(r), by 



h(£,q) 



V2n + 1 \2n 



^ . (2£ + l 



jTTj) J d 2 re- lr ^,(r) 



(B6) 

The Hamiltonian being Gaussian, the height-height cor- 
relation function written in Fourier space is 



(h(£,(i)h(e',q')} = (27r) 2 V<5(q + q') 



AX £ + crq 2 + nq 4 ' 



(B7) 

The height- height correlation function Cjk{r) = 
(hj(r)hk(0)) between membrane j at r and membrane 
k at r = is defined as 



2 k R T 



n-l 



7T 2n 



sm 



2^+1 
2n+T 



dq 



A\ e + aq 2 + nq 4 



ttj sm 



2^+1 
2n + l ' 



irk 



(B8) 



This last equation is thus the generalization of Eq. ( A3 ) 
for n membranes. 

Similar but different results have been obtained in the 
contexts of smectic-A films [43] with a different boundary 
condition for the layer n, and of solid-supported multilay- 



ers [44] using a continuous description of Eq. (B2) valid 



for large n (in the continuum limit d — ¥ 0, n — > oo with 
L = nd held fixed). 
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